Chapter 33 



Semiclassical evolution 



William Rowan Hamilton was born in 1805. At three 
he could read English; by four he began to read Latin, 
Greek and Hebrew, by ten he read Sanskrit, Persian, Ara- 
bic, Chaldee, Syrian and sundry Indian dialects. At age 
seventeen he began to think about optics, and worked out 
his great principle of "Characteristic Function." 
— TurnbuU, Lives of Mathematicians 

(G. Vattay, G. Tanner and P. Cvitanovic) 

SEMICLASSICAL APPROXIMATIONS to quantuHi mcchanics are valid in the regime 
where the de Broglie wavelength A ~ hip of a particle with momentum 
p is much shorter than the length scales across which the potential of the 
system changes significantly. In the short wavelength approximation the particle 
is a point-like object bouncing off potential walls, the same way it does in the 
classical mechanics. The novelty of quantum mechanics is the interference of the 
point-like particle with other versions of itself traveUng along different classical 
trajectories, a feat impossible in classical mechanics. The short wavelength - or remark 33.1 
semiclassical - formalism is developed by formally taking the limit fi ^ m 
quantum mechanics in such a way that quantum quantities go to their classical 
counterparts. 



33.1 Hamilton- Jacobi theory 

We saw in chapter 32 that for a 1-dof particle moving in a slowly varying poten- 
tial, it makes sense to generalize the free particle wave function (32.1) to a wave 
function 

il/iq,t)=Aiq,t)e"^^''''^'^, (33.1) 
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with slowly varying (real) amplitude A{q, t) and rapidly varying (real) phase R{q, t). 

its phase and magnitude. The time evolution of the phase and the magnitude of \p exercise 32.1 

follows from the Schrodinger equation (31.1) 



(33.2) 



Assume A 0, and separate out the real and the imaginary parts. We get two 
equations: The real part governs the time evolution of the phase 



dR 



1 IdR 



dt 2m \ dq 



1 



^ + :^hr +Viq)-—-—A = 0, 



2m A dq^ 



(33.3) 



and the imaginary part the time evolution of the amplitude 



exercise 33.6 
exercise 33.7 



dA I ^ dAdR 1 d^R 
_i_ \ _|_ 

dt m ^ dqi dqi 2m dq^ 



0. 



(33.4) 



In this way a linear PDE for a complex wave function is converted into a set of 
coupled non-linear PDE's for real-valued functions R and A. The coupling term 
in (33.3) is, however, of order fi^ and thus small in the semiclassical limit ^ ^ 0. 

Now we generalize the Wentzel-Kramers-Brillouin (WKB) ansatz for 1-dof 
dynamics to the Van Vleck ansatz in arbitrary dimension: we assume the mag- 
nitude A{q,t) varies slowly compared to the phase R{q,t)lfi, so we drop the fi- 
dependent term. In this approximation the phase R{q, t) and the corresponding 
"momentum field" ^{q, t) can be determined from the amplitude independent 
equation 



exercise 33.8 



OR dR\ ^ 



(33.5) 



In classical mechanics this equation is known as the Hamilton-Jacobi equation. 
We will refer to this step (as well as all leading order in fi approximations to 
follow) as the semiclassical approximation to wave mechanics, and from now on 
work only within this approximation. 



33.1.1 Hamilton's equations 



We now solve the nonlinear partial differential equation (33.5) in a way the 17 
year old Hamilton might have solved it. The main step is the step leading from 
the nonlinear PDE (33.9) to Hamilton's ODEs (33.10). If you already understand 
the Hamilton-Jacobi theory, you can safely skip this section. 
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Figure 33.1: (a) A phase R(q,t) plotted as a 
function of the position q for two infinitesimally 
close times, (b) The phase R(q, t) transported by 
a swarm of "particles"; The Hamilton's equations 
(33.10) construct R{q, t) by transporting ~* 9(0 
and the slope of ^(^o, to), that is po — > pit)- 




(a) 



fast track: 

sect. 33.1.3, p. 639 



The wave equation (31.1) describes how the wave function ij/ evolves with 
time, and if you think of ij/ as an (infinite dimensional) vector, position q plays a 
role of an index. In one spatial dimension the phase R plotted as a function of the 
position q for two different times looks something like figure 33.1 (a): The phase 
R{q, to) deforms smoothly with time into the phase R{q, t) at time t. Hamilton's 
idea was to let a swarm of particles transport R and its slope dR/dq at q at initial 
time t = to to a corresponding R{q, t) and its slope at time t, figure 33.1 (b). For 
notational convenience, define 



dR 

Pi = Piiq,t) := — , i^\,l,...,D. 
oqi 



(33.6) 



We saw earlier that (33.3) reduces in the semiclassical approximation to the Hamilton- 
Jacobi equation (33.5). To make life simple, we shall assume throughout this 
chapter that the Hamilton's function H{q, p) does not depend explicitly on time t, 
i.e., the energy is conserved. 

To start with, we also assume that the function R{q, t) is smooth and well 
defined for every q at the initial time t. This is true for sufficiently short times; 
as we will see later, R develops folds and becomes multi-valued as t progresses. 
Consider now the variation of the function R{q, t) with respect to independent 
infinitesimal variations of the time and space coordinates dt and dq, figure 33.1 (a) 



dR , dR , 
dR = —dt + —dq . 
at dq 



(33.7) 



Dividing through by dt and substituting (33.5) we obtain the total derivative of 
R{q, t) with respect to time along the as yet arbitrary direction q, that is, 



dR 

—{q, q, t) ^ -H{q, p) + q ■ p . 
dt 



(33.8) 



VanVleck - 8dec2010 



ChaosBook.org versionl4.4.I, Apr 21 2013 



CHAPTER 33. SEMICLASSICAL EVOLUTION 



637 



Note that the "momentum" p - dR/dq is a well defined function of q and t. 
In order to integrate R{q, t) with the help of (33.8) we also need to know how 
p = dR/dq changes along q. Varying p with respect to independent infinitesimal 
variations dt and dq and substituting the Hamilton- Jacobi equation (33.5) yields 

dR d^R , d^R , (dH dH dp\ , dp , 
d— - „ „ dt + -r-^dq - - — — h — — ■T-\dt + —dq . 



dq dqdt dq^ \dq dp dq j dq 



Note that H{q,p) depends on q also thi^ough p{q,t) - dR/dq, hence the ^ term 
in the above equation. Dividing again through by dt we get the time derivative of 
dR/dq, that is, 



dH I . dH\ dp 



Piq,q,t) + — = \q-—]^. (33.9) 
dq \ dp j dq 



Time variation of p depends not only on the yet unknown q, but also on the second 
derivatives of R with respect to q with yet unknown time dependence. However, if 
we choose q (which was arbitrary, so far) such that the right hand side of the above 
equation vanishes, we can calculate the function R{q, t) along a specific trajectory 
(^(0, p{t)) given by integrating the ordinary differential equations 

. dH{q,p) . dH{q,p) 

? = ^ , P = 7 (33.10) 

dp dq 

with initial conditions 

dR 

q{to) - q', pito) =p' = —iq', to). (33. 11) 

dq 

We recognize (33.10) as Hamilton's equations of motion of classical mechanics. 
The miracle happens in the step leading from (33.5) to (33.9) - if you missed it, 
you have missed the point. Hamilton derived his equations contemplating optics 
- it took him three more years to realize that all of Newtonian dynamics can be 
profitably recast in this form. 

q is no longer an independent function, and the phase R{q, t) can now be com- 
puted by integrating equation (33.8) along the trajectory {q{t),p{t)) 



R(q,t) = R(q',to)+R(q,t;q',to) 



Riq 



, t; q', fo) = { dT [q{T) ■ p{T) - H{q{T), p(t))] , (33.12) 



with the initial conditions (33.11). In this way the Hamilton- Jacobi partial differ- 
ential equation (33.3) is solved by integrating a set of ordinary differential equa- 
tions, Hamilton's equations. In order to determine R(q, t) for arbitrary q and t we 



section 7.1 



VanVlcck - 8dec2010 



ChaosBook.org vcrsionl4.4.I, Apr 21 2013 



CHAPTER 33. SEMICLASSICAL EVOLUTION 



638 



have to find a q' such that the trajectory starting in {q',p' - dqR{q' , t^)) reaches 
q in time t and then compute R along this trajectory, see figure 33.1 (b). The 
integrand of (33.12) is known as the Lagrangian, 

L{q,q,t) = q ■ p - H{q, p,t) . (33.13) 

A variational principle lurks here, but we shall not make much fuss about it as yet. 

Throughout this chapter we assume that the energy is conserved, and that the 
only time dependence of H{q, p) is through (^(t), p{t)), so the value of R{q, t; q', to) 
does not depend on to, but only on the elapsed time t-to- To simplify notation we 
will set to = and write 

R(q,q',t) = R(q,t;q',0). 



The initial momentum of the particle must coincide with the initial momentum of 
the trajectory connecting q' and q: 

exercise 33.5 

The function R{q, q' , t) is known as Hamilton's principal function. exercise 33.9 

To summarize: Hamilton's achievement was to trade in the Hamilton-Jacobi 
partial differential equation (33.5) describing the evolution of a wave front for a 
finite number of ordinary differential equations of motion, with the initial phase 
R{q, 0) incremented by the integral (33.12) evaluated along the phase space trajec- 
tory {q{T), p{t)). 



33.1.2 Action 

Before proceeding, we note in passing a few facts about Hamiltonian dynamics 
that will be needed for the construction of semiclassical Green's functions. If the 
energy is conserved, the j H{q,p)dT integral in (33.12) is simply Et. The first 
term, or the action 

S{q,q',E)= f dTq{T)-p{T)= {'dq-p (33.15) 

Jo Jq' 

is integrated along a trajectory from q' to q with a fixed energy E. By (33.12) the 
action is a Legendre transform of Hamilton's principal function 

S{q,q',E) = R{q,q',t) + Et. (33.16) 
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The time of flight t along the trajectory connecting q' ^ q with fixed energy E is 
given by 

—S{q,q',E) = t. (33.17) 

The way to think about the formula (33.16) for action is that the time of flight is a 
function of the energy, t = t(q, q', E). The left hand side is explicitly a function of 
E; the right hand side is an implicit function of E through energy dependence of 
the flight time t. 

Going in the opposite direction, the energy of a trajectory E = E{q,q',t) 
connecting q' ^ q with a given time of flight t is given by the derivative of 
Hamilton's principal function 



^R{q,q',t) = -E, (33.18) 
ot 

and the second variations of R and S are related in the standard way of Legendre 
transforms: 



-^Riq,q',t)—Siq,q',E) = -l. (33.19) 

A geometric visualization of what the phase evolution looks like is very helpful 
in understanding the origin of topological indices to be introduced in what fol- 
lows. Given an initial phase R{q,to), the gradient dqR defines a D-dimensional section 33.1.4 
Lagrangian manifold {q, p = dqR(q)) in the full 2d dimensional phase space 
(q, p). The defining property of this manifold is that any contractible loop y in 
it has zero action. 




a fact that follows from the definition of as a gradient, and the Stokes theorem. 
Hamilton's equations of motion preserve this property and map a Lagrangian man- 
ifold into a Lagrangian manifold at a later time, t 

Returning back to the main line of our argument: so far we have determined 
the wave function phase R{q, t). Next we show that the velocity field given by the 
Hamilton's equations together with the continuity equation determines the ampli- 
tude of the wave function. 



33.1.3 Density evolution 

To obtain the full solution of the Schrodinger equation (31.1), we also have to 
integrate (33.4). 

p{q,t) ■-A^ = r^ 
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plays the role of a density. To the leding order in Ji, the gradient of R may be 
interpreted as the semiclassical momentum density 



d dA dR 

ipiqj) {-m—)i//{q,t) = -inA—+p—. 

oq oq oq 



Evaluated along the trajectory (^(0, p{t)), the amplitude equation (33.4) is equiv- 
alent to the continuity equation (16.35) after multiplying (33.4) by 2A, that is 



dp d 

^ + ^(pvO = 0. (33.20) 
ot dqt 



Here, = qi = pi/m denotes a velocity field, which is in turn determined by the 
gradient of R(q, t), or the Lagrangian manifold (q(t), pit) = dqR{q, t)). 



v=-^R{q,t). 
m oq 



As we already know how to solve the Hamilton-Jacobi equation (33.5), we can 
also solve for the density evolution as follows: 

The density p{q) can be visualized as the density of a configuration space 
flow q{t) of a swarm of hypothetical particles; the trajectories q{t) are solutions 
of Hamilton's equations with initial conditions given by (^(0) = q',p{0) - p' - 
d,R{q',0)). 

If we take a small configuration space volume d^q around some point q at time 
t, then the number of particles in it is p{q, t)d^dq. They started initially in a small 
volume d^q' around the point q' of the configuration space. For the moment, we 
assume that there is only one solution, the case of several paths will be considered 
below. The number of particles at time t in the volume is the same as the number 
of particles in the initial volume at ? = 0, 

p{q{t), t)d^q = p{q',0)d^q' , 

see figure 33.2. The ratio of the initial and the final volumes can be expressed as 



p{q{t), t) 



det^ 
oq 



p{q',0). (33.21) 



section 16.2 



As we know how to compute trajectories (^(0, p(t)), we know how to compute 
this Jacobian and, by (33.21), the density p(q(t), t) at time t. 
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Figure 33.2: Density evolution of an initial surface 
(q',p' = dqRiq' ,0) into (q(t), p(t)) surface time t later, 
sketched in 1 dimension. While the number of trajec- 
tories and the phase space Liouville volume are con- 
served, the density of trajectories projected on the q 
coordinate varies; trajectories which started in dq' at 
time zero end up in the interval dq. 



1 



33.1.4 Semiclassical wave function 

Now we have all ingredients to write down the semiclassical wave function at 
time t. Consider first the case when our initial wave function can be written in 
terms of single- valued functions A(q', 0) and Riq', 0). For sufficiently short times, 
R(q, t) will remain a single- valued function of q, and every d^q configuration 
space volume element keeps its orientation. The evolved wave function is in the 
semiclassical approximation then given by 



As the time progresses the Lagrangian manifold dqR{q, t) can develop folds, so 
for longer times the value of the phase R{q, t) is not necessarily unique; in gen- 
eral more than one trajectory will connect points q and q' with different phases 
R{q, q' , t) accumulated along these paths, see figure 33.3. 

We thus expect in general a collection of different trajectories from q' to q 
which we will index by j, with different phase increments Rj{q,q',t). The hy- 
pothetical particles of the density flow at a given configuration space point can 
move with different momenta p = dqRjiq, t). This is not an ambiguity, since in 
the full {q, p) phase space each particle follows its own trajectory with a unique 
momentum. 

Whenever the Lagrangian manifold develops a fold, the density of the phase 
space trajectories in the fold projected on the configuration coordinates diverges. 
As illustrated in figure 33.3, when the Lagrangian manifold develops a fold at 
q = q\; the volume element dq\ in the neighborhood of the folding point is pro- 
portional to ^]dq' instead of dq' . The Jacobian dq' jdq diverges like 1 / ■s[q\'^^'q(7) 
when computed along the trajectory going trough the folding point at ^i. After 
the folding the orientation of the interval dq' has changed when being mapped 
into dq2', in addition the function R, as well as its derivative which defines the 
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Figure 33.3: Folding of the Lagrangian surface 
(q,d,R(q,t)). 



Lagrangian manifold, becomes multi-valued. Distinct trajectories starting from 
different initial points q' can now reach the same final point q2- (That is, the point 
q' may have more than one pre-image.) The projection of a simple fold, or of an 
envelope of a family of phase space trajectories, is called a caustic; this expres- 
sion comes from the Greek word for "capable of burning," evoking the luminous 
patterns that one observes swirling across the bottom of a swimming pool. 

The folding also changes the orientation of the pieces of the Lagrangian man- 
ifold (q, dqR{q, t)) with respect to the initial manifold, so the eigenvalues of the 
Jacobian determinant change sign at each fold crossing. We can keep track of the 
signs by writing the Jacobian determinant as 



det — 

dq 



^~mmj{q,q' ,t) 



det 



dq 



where mj{q, q' , t) counts the number of sign changes of the Jacobian determinant 
on the way from q' to q along the trajectory indexed with j, see figure 33.3. We 
shall refer to the integer mj{q,q',t) as the topological of the trajectory. So in 
general the semiclassical approximation to the wave function is thus a sum over 
possible trajectories that start at any inital q' and end in q in time t 



ilfsc{q,t)= I dq'Y det ^ g.fi,(^,<?',/)/r™,(M',0/2^(y.^0) , (33.22) 
J ^ dq j J 



each contribution weighted by coiTcsponding density, phase increment and the 
topological index. 

That the correct topological index is obtained by simply counting the number 
of eigenvalue sign changes and taking the square root is not obvious - the careful 
argument requires that quantum wave functions evaluated across the folds remain 
single valued. 
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33.2 Semiclassical propagator 

We saw in chapter 3 1 that the evolution of an initial wave function ip(q, 0) is com- 
pletely determined by the propagator (31.10). As K{q,q',t) itself satisfies the 
Schrodinger equation (31.12), we can treat it as a wave function parameterized 
by the configuration point q'. In order to obtain a semiclassical approximation to 
the propagator we follow now the ideas developed in the last section. There is, 
however, one small complication: the initial condition (31.13) demands that the 
propagator at f = is a (5-function at q = q', that is, the amplitude is infinite at 
q' and the phase is not well defined. Our hypothetical cloud of particles is thus 
initially localized at ^ = ^7' with any initial velocity. This is in contrast to the situ- 
ation in the previous section where we assumed that the particles at a given point q 
have well defined velocity (or a discrete set of velocities) given by ^7 = dpH{q, p). 
We will now derive at a semiclassical expression for K{q, q' , t) by considering the 
propagator for short times first, and extrapolating from there to ai^bitrary times t. 

33.2.1 Short time propagator 

For infinitesimally short times 6t away from the singular point ? = we assume 
that it is again possible to write the propagator in terms of a well defined phase 
and amplitude, that is 

K{q, q', 6t) = A(q, q' , ^O^^^^'^'^''*^ . 

As all particles start a.tq = q', R{q, q', 6t) will be of the form (33. 12), that is 

R{q, q, 6t) = pq6t - H{q, p)6t , (33.23) 

with q ~ (q-q')/6t. For Hamiltonians of the form (31.2) we have q = p/m, which 
leads to 

R(q,q',St)^"-^^^^-V{q)dt. 

Here V can be evaluated any place along the trajectory from q to q', for example 
at the midway point V{{q + q')/2). Inserting this into our ansatz for the propagator 
we obtain 

Ksciq, q', 6t) « A{q, q', (JO^sCm )--va/)<sO _ (33 24) 

For infinitesimal times we can neglect the term V{q)dt, so Ksciq,q' ,St) is a d- 
dimensional Gaussian with width cr^ = ifi6t/m. This Gaussian is a finite width 
approximation to the Dirac delta function 

S{z) = lim — ^e-^'/2^' (33.25) 
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if A = (m/lniMt)^^^, with A{q,q' ,6t) fixed by the Dirac delta function normal- 
ization condition. The coiTcctly normalized propagator for infinitesimal times 6t exercise 33.1 
is therefore 

'^^^ ^ [l^r ■ (33-26) 

The short time dynamics of the Lagrangian manifold {q, dqR) which corresponds 
to the quantum propagator can now be deduced from (33.23); one obtains 

dR m 

i.e., is the particles start for short times on a Lagrangian manifold which is a plane 
in phase space, see figure 33.4. Note, that for 6t — > 0, this plane is given by 
the condition q = q' , that is, particles start on a plane parallel to the momentum 
axis. As we have already noted, all particles start at q - q' but with diff"erent 
velocities for ? = 0. The initial surface {q' ,p' = dqR(q',0)) is mapped into the 
surface iq{t), p{t)) some time t later. The slope of the Lagrangian plane for a short 
finite time is given as 

dpi_ _ d^R _ _dp^ _ 
dqj dqjdq'- dqj 5t 

The prefactor (m/6t)^^^ in (33.26) can therefore be interpreted as the determinant 
of the Jacobian of the transformation from final position coordinates q to initial 
momentum coordinates p', that is 

1 / a '\^/^ 

^ I j^^"P \ JR{q,q',dt)/n 



where 



dp[ 
dqj 



t,q' 



d^R{q,q',6t) 
dqjdq'i 



(33.28) 



The subscript • • -1,^/ indicates that the partial derivatives are to be evaluated with 
t, q' fixed. 

The propagator in (33.27) has been obtained for short times. It is, however, 
already more or less in its final form. We only have to evolve our short time 
approximation of the propagator according to (33.22) 

KAq", q',t'+ St) = J]det^ e'^M" ^^M^-i-mM" ^^^<')I^K{q, q'., 6t) , 



J 
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Figure 33.4: Evolution of the semiclassical propaga- 
tor. The configuration which corresponds to the initial 
conditions of the propagator is a Lagrangian manifold 
q = q', that is, a plane parallel to the p axis. The hy- 
pothetical particles are thus initially all placed at q' but 
take on all possible momenta p'. The Jacobian matrix 
C (33.29) relates an initial volume element in momen- 
tum space dp' to a final configuration space volume 
dq. 




and we included here already the possibility that the phase becomes multi-valued, 
that is, that there is more than one path from q' to q". The topological index nij = 
nijiq" ,q' ,t) is the number of singularities in the Jacobian along the trajectory j 
from q' to q". We can write Ksdq" , q' j' + 6t) in closed form using the fact that 
R{q" ,q,t') + R{q,q' ,5t) = R(q",q',t' + St) and the multiplicativity of Jacobian 
determinants, that is 



det 



dq 
dq" 



det 



dq 



dp' 



= det 

g\St dq 



(33.29) 



q',t'+St 



The final form of the semiclassical or Van Vleck propagator, is thus 



det 



dp' 
dq 



1/2 



JRjiq,q',t)/n-imj7i/2 



(33.30) 



This Van Vleck propagator is the essential ingredient of the semiclassical quanti- 
zation to follow. 

The apparent simplicity of the semiclassical propagator is deceptive. The 
wave function is not evolved simply by multiplying by a complex number of mag- 
nitude -y^det dp' /dq and phase R{q, q' , t); the more difficult task in general is to 
find the trajectories connecting q' and ^ in a given time t. 

In addition, we have to treat the approximate propagator (33.30) with some 
care. Unlike the full quantum propagator, which satisfies the group property 
(31.11) exactly, the semiclassical propagator performs this only approximately, 
that is 

Ksciq, q, ti +t2)~ J dq" K^dq, q" , tiW^dq" , q , t\) . (33.31) 

The connection can be made explicit by the stationary phase approximation, sect. 32.2. 
Approximating the integral in (33.31) by integrating only over regions near points 
^" at which the phase is stationary, leads to the stationary phase condition 



dR{q,q",t2) ^ dR{q",q',ti) 



dq" 



dq'.' 



(33.32) 
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Classical trajectories contribute whenever the final momentum for a path from q' 
to q" and the initial momentum for a path from q" to q coincide. Unlike the 
classical evolution of sect. 17.2, the semiclassical evolution is not an evolution by 
linear operator multiplication, but evolution supplemented by a stationary phase 
condition Pf^m = Pm that matches up the classical momenta at each evolution 
step. 



33.2.2 Free particle propagator 

To develop some intuition about the above formalism, consider the case of a free 
particle. For a free particle the potential energy vanishes, the kinetic energy is 
'4q^, and the Hamilton's principal function (33.12) is 



R{q,q',t) = 



/^2 



m{q - q ') 
It 



(33.33) 



The wei ght det ^ from (33.28) can be evaluated explicitly, and the Van Vleck 
propagator is 



D/2 



im(q-q'Yim 



(33.34) 



identical to the short time propagator (33.26), with Viq) = 0. This case is rather 
exceptional: for a free particle the semiclassical propagator turns out to be the 
exact quantum propagator K{q,q',t), as can be checked by substitution in the 
Schrodinger equation (33.2). The Feynman path integral formalism uses this fact remark 33.3 
to construct an exact quantum propagator by integrating the free particle propaga- 
tor (with V{q) treated as constant for short times) along all possible (not necessai^- 
ily classical) paths from q' to q. exercise 33.10 

exercise 33.1 1 
exercise 33.12 



33.3 Semiclassical Green's function 



So far we have derived semiclassical formulas for the time evolution of wave func- 
tions, that is, we obtained approximate solutions to the time dependent Schrodinger 
equation (31.1). Even though we assumed in the calculation a time-independent 
Hamiltonian of the special form (31.2), the derivation would lead to the same final 
result (33.30) were one to consider more complicated or explicitly time dependent 
Hamiltonians. The propagator is thus important when we are interested in finite 
time quantum mechanical effects. For time-independent Hamiltonians, the time 
dependence of the propagator as well as of wave functions is, however, essentially 
given in terms of the energy eigen-spectrum of the system, as in (31.8). It is there- 
fore advantageous to switch from a time representation to an energy representa- 
tion, that is from the propagator (31.10) to the energy dependent Green's function 
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(31.14). A semiclassical approximation of the Green's function Gsc{q,q' ,E) is 
given by the Laplace transform (31.14) of the Van Vleck propagator Ksdq, q' ,t): 

1 r°° 

GAq,q',E) = - dt e'^'"'KAq,q' J) ■ (33.35) 
in Jo 



The expression as it stands is not very useful; in order to evaluate the integral, at 
least to the leading order in h, we need to turn to the method of stationary phase 
again. 



33.3.1 Stationary phase in higher dimensions 

exercise 32.1 

Generalizing the method of sect. 32.2 to d dimensions, consider stationary phase 
points fulfilUng 



d 

— cl)(x) 

dxi 



= V/=1,...J. 

x=xo 



An expansion of the phase up to second order involves now the symmetric matrix 
of second derivatives of 0{x), that is 



Dijixo) - ^^1>(x) 

OXiOXj 



x=xo 



After choosing a suitable coordinate system which diagonalizes D, we can ap- 
proximate the (i-dimensional integral by d 1 -dimensional Fresnel integrals; the 
stationary phase estimate of (32.13) is then 

/ « i^^'/')'"^ Idet D{xoT^^^A{xo) e"<i'(^o)-f , (33.36) 



where the sum runs over all stationary phase points xq of (^{x) and m(xo) counts 
the number of negative eigenvalues of D(xo). exercise 28.2 

exercise 33.2 

The stationary phase approximation is all that is needed for the semiclassical exercise 32.3 
approximation, with the proviso that D in (33.36) has no zero eigenvalues. 



33.3.2 Long trajectories 

When evaluating the integral (33.35) approximately we have to distinguish be- 
tween two types of contributions: those coming from stationary points of the phase 
and those coming from infinitesimally short times. The first type of contributions 



VanVleck - 8dcc2010 



ChaosBook.org versioiiI4.4.1, Apr 21 2013 



CHAPTER 33. SEMICLASSICAL EVOLUTION 



648 



can be obtained by the stationary phase approximation and will be treated in this 
section. The latter originate from the singular behavior of the propagator for f — > 
where the assumption that the amplitude changes slowly compared to the phase 
is not valid. The short time contributions therefore have to be treated separately, 
which we will do in sect. 33.3.3. 

The stationary phase points t* of the integrand in (33.35) ai^e given by the 
condition 



^R{q,q',f) + E = 0. 
ot 



(33.37) 



We recognize this condition as the solution of (33. 18), the time f = f{q, q' , E) in 
which a particle of energy E starting out in q' reaches q. Taking into account the 
second derivative of the phase evaluated at the stationary phase point, 

1 52 

R{q, q', t) + Et = R{q, q', f) + Ef + -{t - ff^R{q, q, f) + ■ ■ ■ 

2 err 

the stationary phase approximation of the integral corresponding to a classical 
trajectory j in the Van Vleck propagator sum (33.30) yields 



Gj{q, q',E) 



1 



detC 



d^Ri 



1/2 



(33.38) 



where nij - mj{q, q' , E) now includes a possible additional phase arising from the 
time stationary phase integration (32.16), and Cj = C j{q,q' ,t*), Rj = Rj{q,q',t*) 
are evaluated at the transit time t* . We re-express the phase in terms of the energy 
dependent action (33.16) 



S{q,q',E) = R{q,q',t*) + Et* , with t* = t*iq,q' , E) , 



(33.39) 



the Legendre transform of Hamilton's principal function. Note that the partial 
derivative of the action (33.39) with respect to qj 

dS{q,q',E) dR{q,q',t*) (dR(q,q',t) \ dt 

— h I h E I . 

dqi dqi \ dt* j dqi 

is equal to 



dS{q,q',E) dR{q,q',f) 



dqi 



dqi 



(33.40) 



due to the stationary phase condition (33.37), so the definition of momentum as a 
partial derivative with respect to q remains unaltered by the Legendre transform 
from time to energy domain. exercise 33.13 
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Next we will simplify the amplitude term in (33.38) and rewrite it as an ex- 
plicit function of the energy. Consider the [(D + l)x(D + 1)] matrix 



dqdE dE^ 




(33.41) 



where S = S{q,q' , E) and we used (33.14-33.17) here to obtain the left hand side 
of (33.41). The minus signs follow from the definition of (33.15), which implies 
that that S (q, q',E) = -S (q',q, E). Note that D is nothing but the Jacobian matrix 
of the coordinate transformation {q,E) (p'j) for fixed q'. We can therefore 
use the multiplication rules of determinants of Jacobians, which are just ratios of 
volume elements, to obtain 



We use here the notation (det .)^y , for a Jacobian determinant with partial deriva- 
tives evaluated at t, q' fixed, and likewise for other subscripts. Using the relation 
(33. 19) which relates the term ^ to d^R we can write the determinant of D as 
a product of the Van Vleck determinant (33.28) and the amplitude factor arising 
from the stationary phase approximation. The amplitude in (33.38) can thus be 
interpreted as the determinant of a Jacobian of a coordinate transformation which 
includes time and energy as independent coordinates. This causes the increase in 
the dimensionality of the matrix D relative to the Van Vleck determinant (33.28). 

We can now write down the semiclassical approximation of the contribution 
of the jth trajectory to the Green's function (33.38) in explicitly energy dependent 
form: 



However, this is still not the most convenient form of the Green's function. 

The trajectory contributing to Gj{q,q',E) is constrained to a given energy 
E, and will therefore be on a phase space manifold of constant energy, that is 
H(q, p) = E. Writing this condition as a partial differential equation for S (q, q',E), 
that is 



H{q,—)^E, 
oq 
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one obtains 



^H(q,p) =0 = 
dq. 

^H{q',p') =0 = 
oqi 



dH dpj 
dpj dq'. 



1j 



dqjdq'. 



dqidq'. 



(33.43) 



that is the sub-matrix d^S I dqidq' j has (left- and right-) eigenvectors corresponding 
to an eigenvalue 0. Rotate the local coordinate system at the either end of the 
trajectory 

{q\,q2,qj„ - ■ ■ ,qd) {q\\,q±i,q±2, - ■ ■ ,q±(D-i)) 

so that one axis points along the trajectory and all others aie perpendicular to it 

iqi,q2, 43, ■■■,qd) ^ (q, 0, 0, • • • , 0) . 

With such local coordinate systems at both ends, with the longitudinal coordinate 
axis q\\ pointing along the velocity vector of magnitude q, the stability matrix of 
S{q, q', E) has a column and a row of zeros as (33.43) takes form 



d^S 



d^S 



dq\\dq'. dqidq 



jq' = 0. 



The initial and final velocities are non-vanishing except for points \q\ = 0. These 
are the turning points (where all energy is potential), and we assume that neither q 
nor q' is a turning point (in our application - periodic orbits - we can always chose 
q = q' not a turning point). In the local coordinate system with one axis along 
the trajectory and all other perpendicular to it the determinant of (33.41) is of the 
form 



detD(^,^',£) = (-1) 



D+l 





det 





d- s 



dq\\dE 



d-s 

dEM. 



(33.44) 



The corner entries can be evaluated using (33.17) 



d^S 



—t-- 



d^S 



1 



dq\\dE dq\\ q dEdq'. q' 



As the q\\ axis points along the velocity direction, velocities q, q' are by construc- 
tion almost always positive non-vanishing numbers. In this way the determinant 
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of the [(D + l)x(D + 1)] dimensional matrix D{q,q' ,E) can be reduced to the 
determinant of a [(D - l)x(D - 1)] dimensional transverse matrix Dj^{q,q',E) 



det D{q, q',E) ^ — det D_L_{q, q' , E) 
qq' 

n r ' d^S{q,q',E) 
Dj_{q,q,E)ik ^ 



dq±idq' 



(33.45) 



±k 



Putting everything together we obtain the jth trajectory contribution to the semi- 
classical Green's function exercise 33.15 



Gj{q,q',E) = 



1/2 



detD^ 



1/2 



(33.46) 



where the topological index mj = mj(q, q', E) now counts the number of changes 
of sign of detD^ along the trajectory j which connects q' to q at energy E. 



The endpoint velocities q, q' also depend on {q, q' ,E) and the trajectory / 



33.3.3 Short trajectories 

The stationary phase method cannot be used when t* is small, both because we 
cannot extend the integration in (32.16) to -oo, and because the amplitude of 
K{q,q',t) is divergent. In this case we have to evaluate the integral involving the 
short time form of the exact quantum mechanical propagator (33.26) 

Goiq,q',E) = i j^"^ . (33.47) 



By introducing a dimensionless variable t - t^2m(E - V(q))/rn\q - q'\, the 
integral can be rewritten as 



where Soiq, q' ,E) = ^J2m{E - V)\q - q'\ i& the short distance form of the action. 
Using the integral representation of the Hankel function of first kind 



" Jo 



^ife(r+l/r)^-y-l^^ 
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we can write the short distance form of the Green's function as 



im I ^J2m{E-V)\ 
Go(q,q',E)^-—{-^ HU{So{q,q' ,E)in) . (33.48) 



Iffi \ 2nn\q - q'\ 

Hankel functions are stabdai^d, and their the short wavelength asymptotics is de- 
scribed in standard reference books. The short distance Green's function approx- 
imation is vahd when So{q, q',E) < fi. 



Resume 



The aim of the semiclassical or short-wavelength methods is to approximate a 
solution of the Schrodinger equation with a semiclassical wave function 



j 



iRj{q,t)/n 



accurate to the leading order in H. Here the sum is over all classical trajectories 
that connect the initial point q' to the final point q in time t. "Semi-" refers to h, 
the quantum unit of phase in the exponent. The quantum mechanics enters only 
through this atomic scale, in units of which the variation of the phase across the 
classical potential is assumed to be large, "-classical" refers to the rest - both the 
amplitudes Aj{q, t) and the phases Rj{q, t) - which are determined by the classical 
Hamilton- Jacobi equations. 

In the semiclassical approximation the quantum time evolution operator is 
given by the semiclassical propagator 



KsM,q',t) = 



1 



det 



dp' 
dq 



1/2 



where the topological index mj{q, q' , t) counts the number of the direction reversal 
along the j'th classical trajectory that connects q' ^ qin time t. Until very recently 
it was not possible to resolve quantum evolution on quantum time scales (such as 
one revolution of electron around a nucleus) - physical measurements are almost 
always done at time scales asymptotically large compared to the intrinsic quantum 
time scale. Formally this information is extracted by means of a Laplace trans- 
form of the propagator which yields the energy dependent semiclassical Green 's 
function 



Gsc{q,q',E) = GQ{q,q' ,E) + ^G j{q,q' ,E) 

j 

1 1 f>n' 

Gjiq,q',E) 



ifi{2nifi) 



(D-i) 



qq' dq_L 



(33.49) 
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where GQ{q,q' ,E) is the contribution of short trajectories with So(q,q',E) < fi, 
while the sum is over the contributions of long trajectories (33.46) going from q' 
to q with fixed energy E, with S j{q, q',E) :s> fi. 

Commentary 

Remark 33.1 Limit — > 0. The semiclassical limit "S — > 0" discussed in sect. 33.1 
is a shorthand notation for the limit in which typical quantities like the actions R or S in 
semiclassical expressions for the propagator or the Green's function become large com- 
pared to H. In the world that we live in the quantity ?z is a fixed physical constant whose 
value [33.8] is 1.054571596(82) IQ-^"* Js. 

Remark 33.2 Madelung's fluid dynamics. Already Schrodinger [33.3] noted that 



plays the role of a density, and that the gradient of R may be interpreted as a local semi- 
classical momentum, as the momentum density is 



oq oq oq 



A very different interpretation of (33.3-33.4) has been given by Madelung [33.2], and 
then built upon by Bohm [33.6] and others [33.3, 33.7]. Keeping the fi dependent term 
in (33.3), the ordinary differential equations driving the flow (33.10) have to be altered; 
if the Hamiltonian can be written as kinetic plus potential term V{q) as in (3 1 .2), the tP' 
term modifies the p equation of motion as 



-^iViq) + Q(q,t)) , (33.50) 
oqi 



where, for the example at hand, 



1 

Q{q,f)^— pT^^ (33.51) 

2m yfp oq^ 



interpreted by Bohm [33.6] as the "quantum potential." Madelung observed that Hamil- 
ton's equation for the momentum (33.50) can be rewritten as 



^+ — v, = o-y. (33.52) 

at \ oqj m oqi mp oqj 
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where cr, ,■ — 



ti-p d-lnp • 
4m dqjdqj 



is the "pressure" stress tensor, v. 



',■ = Pi/m, and p - as de- 



fined [33.3] in sect. 33.1.3. We recall that the Eulerian ;^ + is the ordinary deriva- 
tive of Lagrangian mechanics, that is ^ . For comparison, the Euler equation for classical 
hydrodynamics is 



where pdij is the pressure tensor. 

The classical dynamics corresponding to quantum evolution is thus that of an "hypo- 
thetical fluid" experiencing h and p dependent stresses. The "hydrodynamic" interpreta- 
tion of quantum mechanics has, however, not been very fruitful in practice. 

Remark 33.3 Path integrals. The semiclassical propagator (33.30) can also be derived 
from Feynman's path integral formalism. Dirac was the first to discover that in the short- 
time limit the quantum propagator (33.34) is exact. Feynman noted in 1946 that one can 
construct the exact propagator of the quantum Schrodinger equation by formally summing 
over all possible (and emphatically not classical) paths from q' io q . 

Gutzwiller started from the path integral to rederive Van Vleck's semiclassical ex- 
pression for the propagator; Van Vleck's original derivation is very much in the spirit 
of what has presented in this chapter. He did, however, not consider the possibility of 
the formation of caustics or folds of Lagrangian manifolds and thus did not include the 
topological phases in his semiclassical expression for the propagator Some 40 years later 
Gutzwiller [34.4] added the topological indices when deriving the semiclassical propaga- 
tor from Feynman's path integral by stationary phase conditions. 

Remark 33.4 Applications of the semiclassical Green's function. The semiclassical 
Green's function is the starting point of the semiclassical approximation in many appli- 
cations. The generic semiclassical strategy is to express physical quantities (for example 
scattering amplitudes and cross section in scattering theory, oscillator strength in spec- 
troscopy, and conductance in mesoscopic physics) in terms of the exact Green's function 
and then replace it with the semiclassical formula. 

Remark 33.5 The quasiclassical approximation The quasiclassical approximation 
was introduced by Maslov [32.7]. The term 'quasiclassical' is more appropriate than 
semiclassical since the Maslov type description leads to a pure classical evolution operator 
in a natural way. Following mostly ref. [32.8], we give a summary of the quasiclassical 
approximation, which was worked out by Maslov [32.7] in this form. One additional 
advantage of this description is that the wave function evolves along one single classical 
trajectory and we do not have to compute sums over increasing numbers of classical 
trajectories as in computations involving Van Vleck formula [34.27]. 
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Exercises 



33.1. Dirac delta function, Gaussian representation. 

Consider the Gaussian distribution function 



33.6. 1 -dimensional harmonic oscillator. 



Take a 1 - 



SAz) 



1 



Show that in cr —> Hmit this is the Dirac delta function 



Jm 



dx 6{x) = 1 if e A1 , zero otherwise . 



33.2. Stationary phase approximation in higher dimensions. 

All semiclassical approximations are based on saddle 
point evaluations of integrals of type 



dimensional harmonic oscillator U(q) = jkq~. Take a 
WKB wave function of form A{q, t) - a(t) and R(q, t) — 
r(t) + b(t)q + c{t)q^, where r(f), a{t), b{t) and c(f) are 
time dependent coefficients. Derive ordinary differen- 
tial equations by using (33.3) and (33.4) and solve them, 
(continuation: exercise 33.9.) 

33.7. 1 -dimensional Unear potential. Take a 1 -dimensional 
linear potential U (q) = -Fq. Take a WKB wave func- 
tion of form A{q, t) = a{t) and R{q, t) - r{t) + b{t)q + 
c(t)q^, where r(f), a(t), b(t) and c(t) are time dependent 
coefficients. Derive and solve the ordinary differential 
equations from (33.3) and (33.4). 



xA(x)e 



i(s>(x)/n 



(33.53) 



for small values of h. Obtain the stationary phase esti- 
mate 



/ ~ ^ A{x„)e 



,o(x„)//,^27n^ 



VdetD2(D(x„) 

where D^<l)(x„) denotes the second derivative matrix. 

33.3. Schrodinger equation in the Madelung form. Ver- 
ify the decomposition of Schrodinger equation into real 
and imaginary parts, eqs. (33.3) and (33.4). 



33.4. Transport equations. 

function in the asymptotic form 



Write the wave- 



il,(qJ)^ei''^'''^^i''Y,m"A„(x,t). 



H>0 



Derive the transport equations for the A„ by substituting 
this into the Schrodinger equation and then collecting 
terms by orders of fi. Note that equation for A„ only 
requires knowledge of A„_i and R. 



33.8. D-dimensional quadratic potentials. Generalize the 
above method to general D-dimensional quadratic po- 
tentials. 

33.9. Time evolution of R. (continuation of exer- 
cise 33.6) Calculate the time evolution of R{q,0) - 
a + bq + cq^ for a 1 -dimensional harmonic oscillator us- 
ing (33.12) and (33.14). 

33.10. D-dimensional free particle propagator. Verify the 
results in sect. 33.2.2; show explicitly that (33.34), the 
semiclassical Van Vleck propagator in D dimensions, 
solves the Schrodinger's equation. 

33.11. Propagator, charged particle in constant magnetic 

field. Calculate the semiclassical propagator for a 
charged particle in constant magnetic field in 3 dimen- 
sions. Verify that the semiclassical expression coincides 
with the exact solution. 

33.12. 1-dimensional harmonic oscillator propagator. 

Calculate the semiclassical propagator for a 1- 
dimensional harmonic oscillator and verify that it is 
identical to the exact quantum propagator. 

33.13. Free particle action. Calculate the energy dependent 
action for a free particle, a charged particle in a constant 
magnetic field and for the harmonic oscillator 



33.5. Easy examples of the Hamilton's principal function. 33.14. 
Calculate R{q, q', t) for 

a) a D-dimensional free particle 

b) a 3 -dimensional particle in constant magnetic field 

c) a 1 -dimensional harmonic oscillator. 



Zero length orbits. Derive the classical 

trace (18.1) rigorously and either add the f — » 0+ zero 
length contribution to the trace formula, or show that it 
vanishes. Send us a reprint of Phys. Rev. Lett, with the 
correct derivation. 



(continuation: exercise 33.13.) 



33.15. Free particle semiclassical Green's functions. Cal- 
culate the semiclassical Green's functions for the sys- 
tems of exercise 33.13. 
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